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ABSTRACT 

We present a new more accurate approach to the composite spectra construction 
based on stacking spectra with similar slopes a\ within the wavelength range redward 
of Lya emission line, which allows to reduce a noise. With the help of this technique 
a detailed study of the Hi Lya-forest region (Arcst ~ 1050 — 1200 A) of the own 
sample of 3,439 medium-resolution quasar spectra from SDSS DR7 was performed. 
More than 14 lines were found within it, three of which were found in previous studies 
of quasar composite spectra from SDSS and some others were found in composite 
spectra from space-based telescopes or high-resolution spectra of individual quasars 
from ground-based telescopes. The parameters of these lines were calculated. It was 
shown that the continuum level within the Lya-forest region cannot be considered 
as a power-low with the same slope as in the region redward of Lya emission line. 
Any dependence of the slope a a on luminosity in SDSS u, g, r and i bands as far as 
on luminosity in 1450 — 1470 A band was not found. The proposed approach can be 
applied for generation of new templates for more precise quasar redshift measurements, 
theoretical determination of K-correction and color-indexes, as far as for determination 
of continuum and mean transmission in Lya forest studies. It was shown that the 
uncertainties in the mean transmission caused by using the composite spectra made 
with common approach can constitute about 20%. 
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1 INTRODUCTION 

The Lya forest in the spectra of distant quasars traces 
the thermal and radiative history of the Universe, as well as 
the evolution of underlying matter distribution over a wide 
range of scales and redshifts. It is possible due to relation 
of the Lya opacity of the intergalactic neutral hydrogen H I 
to its density and other physical parameters. As a measure 
of opacity the value F named the transmission is used; fluc- 
tuations of it, 5f, provide an invaluable information about 
the density fluctuati ons on the smal l est scales available to- 
day for observations jKim et al.ll2004 ■ McDonald et abllioO^ : 



iD'bdorico et al.ll2006l : ICappetta et al.ll20ld ) 

Two important steps in these studies, after the row 
spectral data reduction and before calculation of the 
two-point statistics (correlation function and power spec- 



* E-mail: g.ivashchenkoOgmail.com 
t E-mail: olka@astro.franko.lviv.ua 
X E-mail: el.torbaniuk@gmail.com 



trum) of transmission fluctuations are (i) determination of 
the redshift-dependent mean transmission F level (around 
which one measures the fluctuations) and (ii) determina- 
tion of the intrinsic quasar spectrum (transmission of which 
one measures). For today several tens of high-resolution 
(R ~ 100, 000) quasar spectra and several thousands of 
medium-resolution {R ~ 2, 000) ones are available for Lya- 
forest study. The flrst ones ar e obtained with sp ectrographs 
like UVES on VLT (see e.g. iKim et all !{20o4 )) and allow 
a good spatial resolution for two-point statistics, but using 
them along one cannot be sure about the reliability of the 
mean transmission level for a given redshift due to small 
number of spectra. Moreover in this case the quasar contin- 
uum is fitted manually with the spline polynomials without 
any underlying ideas about the emission lines' shape etc., 
that can result in additional uncertainties. It is worth to 
note, that the term continuum in this case usually means 
the whole intrinsic quasar spectrum including emission lines. 
We follow this definition and refer to the continuum without 
emission lines as the true continuum. The medium-resolution 
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spectra are obtained mainly in the scope of the Sloan Digi- 
tal Sky Survey (SDSS) project ( Abazaiian et al.ji2009 '). They 
can serve well for the mean transmission determination, but 
they do not allow the manual direct fitting of the continuum 
within the Lya-forest part because it is almost impossible 
to find unabsorbed regions in such spectra. 

For calculations of the mean transmission and contin- 
uum level in the case of medium resolution spectra the so- 
called composite (mean) spectra are used. The possibility 
of their utilization is related to the fact that quasar spec- 
tra are remarkably similar from one object to another. Such 
ultraviolet-optical spectra can be characterized by a feature- 
less continuum and a series of mostly broad emission line 
features; and due to their high signal-to-noise ratio these 
spectra also reveal weak emission features that are rarely 
detectable in individual quasar spectra. The composite spec- 
tra were compiled for a wide set of qua sar samples, e. g. th e 
Large Bright Quasar Survey (LBQS) llFrancis et al.l 199 j) 



the F irst Bright Qu asar Survey (FBQS) fBrothcrt on et al 



200ir) . the SDSS (IVanden Berk et al.1 [2001; F ieri et al 



2010f ). quasar spectra from the Hubble S pace Telescope 



(HST) (|Zheng et al.ll 19971 : iTelfer et al.ll20oj) and Far Ult ra- 
violet Spectroscopic Explorer (FUSE) (|Scott et al1l2004l '). 

The part of UV-optical quasar spectrum redward of 
the Lya (1215.67 A) emission line (free of hydrogen Lyman- 
series forests of abso rption lines) is stud i ed qu ite well. As 
it was shown e. g. by IVanden Berk et al.l (|200ll ) this region 
up to ~ 5000A is described well with a smooth power-low 
shap e ~ A"^ of the true con tinuum and broad emission 
lines. IVanden Berk et all (|200j) reported on a\ to be —1.54 
within 1350-4230 A, that agrees w ell with the values —1.54 
and —1.68 which were obtained by i Brotherton et al.1 (|200ll ') 
within 1450-5050 A for spectra from LBQS and FBQS cor- 
respo ndingly, and also with the results of ICarballo et al.l 
1 199911 . who obtained the values -1.34 ± 0.15 and -2.11 ± 
0.16 within 1300-3000 A and 3000-4500 A ranges corre- 
spondingly for a sample of radio-loud quasars. Throughout 
this paper the wavelengths mean rest-frame ones, unless oth- 
erwise is specified. 

The part of the spectrum blueward of Lya emis- 
sion line, or strictly speaking the part of spectrum be- 
tween dublet of Ly/? (1025.72 A) - O vi (1033.83 A) and 
Lya, which is used for the Lya forest region studies, is 
much less better investigated due to the presence of the 
Lya forest itself. There are two techniques used for con- 
tinuum level determination within this region. The first 
one lies in the direct fitting, usually by polynomials, of 
it within the Lya forest region in high-resolution spectra 
(iMcDonald et al.ll2000l: iD'Odorico et al.ll2006l : [Becker et al] 
I2OO7I : iFaucher-Giguere et al.ll2008l ). For medium-resolution 
spectra, where it is difficult to select unabsorbed regions, 
another technique is applied: it is based on extrapolation 
of the true continuum level from the red side of the Lya 
emission line, that is free from H I absorption with addi- 
tion of several emission lines or without it and determina- 
tion of continuum level directly in the Lya region model- 
ing the spectrum with the several-parametric model includ- 
ing the mean transmission and paramete r s of emission lines 
llPress et all Il993l: iBernardi et al.l l2003l: iDesiacaues et al.l 



spectra these authors worked with a whole sample of indi- 
vidual spectra and applying the principal component analy- 
sis found the mean continuum level, mean transmission and 
mean generalized calibration vector, a function defined by 
them to include calibration errors and mean absorption by 
metal lines with A > 1300A. and also their fluctuations. 



Several-parametric modelling of the continuum level 
within the Lya forest region suffers from the degeneration 
between the mean transmission and parameters of the true 
continuum, that is why extrapolation of the true continuum 
from the red part of spectrum seems to be much reliable. 
But in practice it is only an approximation of the true con- 
tirmum because of the fact that the peak of the well-known 
big blue bump in quasar spectra is located within 1000- 
1300 A region. This is confirmed by the values of a\ ob- 
tained from HST and FU SE sp ectra, -1.01 ± 0.05 within 
I1997I ). -0.24 ± 0.12 within 500- 
2003), -l-44l^ ?^ within 630-1155 A 



1050-20 20 A llZheng et al, 
1200 A (iTelfer et al.! ' 



I2OO7I : iPolinovskvi fc MalvginI I2OI0I : ISongailal I2OO4I ). More 
complica ted, but similar t o the last one, method was pro- 
posed bv lMcDonald et"aLl l|200(j) . Instead of using composite 



— — ■ ■ — ■/ ' — u.^o 

l|Scott et al.ll2004l ). when comparing with the same values 
for longer wavelengths mentioned above. Hence the several- 
parametric modelling is left to be the most reliable method. 
But except the degeneration of parameters its usage en- 
counters another problem, the choice of the model for the 
continuum level. Usually a superposition of the power-low 
and three Gaussians is considered. These Gaussians present 
three broad emission features, two around 1073 A and 1123 A 
an d Lya fine. The values of 1073 A and 1 123 A were found 
by Bernardi et al] l|2003l ) and used by IDesiacaues et alj 
l|2007l ). But there are also other estimations of these val- 
ues. E. g. IVanden Berk et alj (|200ll ) measured them to be 
1065.10 A and 1117.26 A from the SDSS quasar composite 
spectrum and referred to them as An [Xiab = 1066.66 A) 
and Felll multiplet; they also found one more emission line. 
Cm* {Xiab = 1175.67 A) with the wavelength 1175.35 A, 
within L yg-forest reg i on. T he same emission features were 
found bv lTvtler et al] (|2004l ) to be 1070.95 A, 1123.13 A and 
1175.88 A correspondingly. More close look to these emission 
features during the traditional modelling mentioned above 
and performed by us showed that the values of these cen- 
tral wavelengths vary significantly from spectrum to spec- 
trum and the 'line' profiles in some cases are far enough 
from symmetric Gaussians. These conclusions together with 
the results on the s pace-based c omposi te UV spectra o f 
the nearest quasars l|Telfer et all [20021 : IScott et al] |2004| ) 
and ground-based individual optical spectra of some quasars 
(|Binette fc Krongold|[2008[ : iLeighlv et a"i]|2007l ) led us to the 
idea of the improvement of the model. One of the main 
sources of uncertainties during the generation of the com- 
posite spectrum lies in a fact that despite of their similarities 
in general the quasar spectra differ in the slopes of the true 
continuum level in the red part of UV-optical bump. Thus 
we tried to generate the composite spectra from the samples 
of spectra with the similar ax- 

The data used in this work, the sample selection and the 
composite spectra generation processes are described in Sec- 
tion[2l The methods of lines identification and calculation of 
their parameters are presented in Section O The results and 
their discussion are presented in Section[4l Finally Section[5] 
contains the discussion of obtained results. 
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Figure 1. Redshift distribution of all objects from the HW sam- 
ple with 2.3 < z < 4.6 and redshift determination confidence level 
> 0.9 (dash-dot), all visually selected (dotted) for the present 
study and those with the rms of the normalization constant A 
less than 15% (solid). 



2 DATA 

2.1 The SDSS DR7 quasar sample 

Our sample is taken from the public available release 
of the sky-residual subtracted spectra for the Sloan DiRi- 
tal S ky Survey (SDSS) Legacy Release 7 (|Abazaiian et al.l 
I2OO9I ). T his release ('WH samp le' hereafter), which is de- 
scribed in lWild fc HewettI (120101'). contains a total of 106,006 
spectra, generated using the lWild fc HewettI (|2005|) scheme 
from the spectra of the objects from the Schneider et al.l 
(|2010l ) quasar catalogue. This scheme includes a significantly 
improved technique of the OH sky lines subtraction. The 
strong OH sky emission lines extend over almost half of the 
wavelength range (> 6700 A) of the SDSS spectra and is 
not subtracted optimally enough with the previous SDSS 
pipelines. 

The redshifts in SDSS are measured mainly with the 
help of two techniques — emission line measurements and 
cross-correlation, or one of them. The errors of redshift (due 
to intrinsic line shifts etc.) could affect the results of our 
study hence we tried to minimize their influe nce using the 
impro ved redshifts of the objects from the .Schneider et al.l 
quasar c atalog ue, which was generated using the 
iHewett fc Wild l|2010h scheme, inst ead of redshifts in th e 
headers of each fits-file. According to lHewett fc WildI (|2010h . 
their redshifts possess systematic biases, which are a factor 
of 20 lower compared to the SDSS redshift values. 



2.2 Spectra selection 

All the quasar spectra (15,154 objects) from the HW 
sample with redshifts within the range 2.3—4.6 and redshift 
determination confidence level > 0.9 were firstly selected. 
Because the SDSS is an automatic survey there is a possi- 
bility of pollution of the sample with 'wrong' objects whose 
general photometry or spectroscopy properties required by 
the automatic selection pipelines are similar to the ones of 
'real' objects. Thus a visual examination was carried out. 
Except the non-quasar spectra the spectra of quasars with 
low signal-to-noise ratio, broad absorption lines (BAL) and 
Damped Lya (DLA) systems were also excluded during this 
examination. 

The usual formal definition of BAL quasar is a quasar 
with a positive balnicity index (a measure of the equivalent 



width of the C iv absorption) jWevmann et al.llioOll 'l . Their 
spectra show broad (~ 2000 — 20000 km/s, sometimes up to 
66,000 km/s), blueshifted absorption troughs corresponding 
to highly-ionized rest-frame-UV resonance transitions such 
as C IV, Si IV, N V, O vi, more ra rely Mg 11 and Al i ll, 
and even more rarely Fe 11 (see e. g. iMenou et al.l HoOl')). 
The similarity of the continuum and line emission of BAL 
and non-BAL quasars motivates the hypothesis that BAL 
quasars are not intrinsically different from other quasars and 
the presence of BAL features in the spectra of only ~ 10% of 
optically selected quasars is usually explained by a difference 
in viewing angle if the subrelativistic outfiow at the origin of 
BAL features is non isotropic. The popular notion is that the 
outfiow is preferentially located in the plane of the disc sur- 
rounding the supermassive black hole. The presence of the 
BAL quasars in the sample could introduce additional errors 
into an estimation of the quasar mean continuum level. 

The DLA systems are a class of quasar absorbers se- 
lected for the presence of H I column densities > 2-10^'' cm~^ 
(see e.g.jWolfc et al. (2005) for review). They are identified 
as absorption features with the rest equivalent width exceed- 
ing 5 A. The nature of these systems is still not understood, 
however due to their high density they are usually related 
to the galaxy formation and therefore could not be used 
as representatives of the linear perturbations in the neutral 
intergalactic medium. 

After the visual examination the sample contains 
4,779 spectra. Redshift distributions of these objects and all 
the objects initially selected from HW sample are presented 
in Figme\T\(dash-dot and dotted correspondingly). 



2.3 Generating composite spectra 

Generating composite spectra of quasars includes three 
steps: (i) normalization of each spectrum, (ii) setting each 
spectrum to the rest frame and (iii) calculation of the mean 
spectrum. Before we smoothed all spectr a with a simple 
movin g average over 3 points. Following iMcDonald et al.l 
(|2006l ) we removed the following wavelength regions from 
our analysis because of calibration problems due to strong 
sky lines: 5575A < A < 5583A, 5888A < A < 5895A, 
6296A < A < 6308A, and 6862A < A < 687lA, 
wh ere the second r e gion i s lA wider than that presented 
in iMcDonald et al. I l|2006l) . because we added also Na I 
(5894.6 A) interstellar line to it. Normalization of each spec- 
trum is needed to be applied due to different apparent flux 
densities. Taking into accoun t the s imil arity of quasa r spec - 
tra and following iPress et al. I l|l993t ) and lZheng et~aLl (ll997D 
we normalize each spectrum on the (arithmetic) mean flux 
in all pixels within the rest wavelength range 1450-1470 A. 
This range lies blueward from the C iv emission line and 
is usually considered to be free of obvious emission and ab- 
sorption. The last claim appears to be not precise because of 
the weak emission line, which can be seen in composite spec- 
tra, but we neglect this fact, because intensity of this line is 
compared with noise level in individual spectra. For further 
study to reduce possible uncertainties we used only spectra 
with the rms of the normalization constant A less than 15%; 
the number of them is 3,493 and its redshift distribution is 
shown in Figure \I\(solid line). 

Considering the true continuum of the quasar spectra 
redward of Lya emission line to be a power- low ~ A"-'' , 
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Figure 2. The distribution of tlie true continuum slopes rcdward 
of hya emission line. 
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Figure 3. Composite spectra of 16 subsamples witliin 1025- 
1650 A. The highest spectrum has the steeper true continuum 
redward of Ly« emission line. 

we calculated its slope ax for each individual spectrum 
witliin the range between Lya and Civ(1549.0 A) emission 
lines. For this purpose we selected the following wavelength 
ranges from the composite spectrum compiled from all the 
spectra: 1278-1286 AA, 1320-1326 AA, 1345-1360 AA and 
1340— 1480 AA. The distribution of the obtained slopes is 
shown in Figure [51 

Using the obtained values of ax , we selected 16 subsam- 
ples each of 200 spectra with ax closest to —0.7 — k ■ 0.1, 
k — 0..15. Then dividing each j-th selected spectrum onto its 
normalization constant A-' , rebinning them with AArest ~ 
2 A and stacking the spectra into the rest frame, we ob- 
tained the mean arithmetic composite spectra. The disper- 
sion of each pixel of the composite spectrum was cal- 
culated from the noises ai of pixels of individual spectra 
constituting it as = f,"^- These spectra within 1025- 

i 

1650 A and enlarged part within 1025-1200 A are shown in 
Figures [3] and [4] correspondingly. The dashed lines indicate 
the laboratory wavelengths of the indicated lines identified 
by [Vandcn Berk et al.. (,2001, ) in this range. The slopes of 
these composite spectra, calculated as mentioned above, 
with the errors of approximation, and also the mean red- 
shift and the mean absolute magnitude in SDSS i-band Mi 
of each subsample with their rms are presented in Table [1] 
The slopes of the true continuum in this case were cal- 
culated in the way mentioned above, but the parts of the 
spectra which are the most free from emission lines were se- 
lected manually in each spectrum. The part of spectra used 
for calculation of the true continuum slope is shown in Fig- 
ure^{short-dashed line, the highest spectra are the steepest 
ones). The three spectra with the slopes —0.84, —1.49 and 
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Figure 4. The same as in Fig. E] within 1025-1200 A. 
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Figure 5. The parts of the composite spectra used for true con- 
tinuum approximation. See explanation in the text. 

— 2.12 are shown with the solid lines together with the fitted 
power- low true continua {dashed lines). The regions which 
were used for fitting the true continuum slightly vary from 
spectrum to spectrum and not all four parts used in calcu- 
lation of the slopes of the individual spectra we taken into 
account in the case of composite spectra. This is obviously 
seen on an example of the region with the shortest wave- 
lengths: with the steepening of the spectrum the centre of 
it moves from ~ 1280 A to ~ 1290 A. It means that approx- 
imation of the true continuum using the same regions as it 
was done for individual spectra is only an approximation, 
but this approximation is sufficient for rough separation of 
spectra for composites compilation. 

The values of absolute magnitudes in i-band. Mi, for 
each individual object were taken from the Sloan D igital Sky 
Survey Quasar Catalog V (jSchneider et al.ll2010l ). where it 
is calculated within the frame of the flat ACDM cosmolog- 
ical model with Ho = 70 km/s/Mpc, Q.u = 0.3, Qa = 0.7, 
and correcting for G alactic extinction (using the maps of 
ISchlegel et alT (|l998h ) and assuming that the quasar spec- 
tral energy distribution in the UV-optical can be represented 
by a power law with ax = —1.5. As one can see, the objects 
in all subsamples have similar mean redshift and mean ab- 
solute magnitude. It means that any possible differences in 
the spectral properties within the wavelength region under 
consideration, which can be obtained in course of the present 
study, cannot be explained by the evolution effects or related 
to quasar luminosity, at least in i-band. 

It is worth to note that the absolute magnitude in i- 
band usually used as a value characterizing the luminos- 
ity of objects is not a good characteristic in case of quasar 
UV-bump because its mean wavelength is 7481 A, which is 
far from the peak of the big blue bump, that is why we 
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Figure 6. Slope ax — absolute magnitude M diagrams for u, g, r, i-bands. 



Table 1. Characteristics of 16 subsamples, the mean redshift 
z and the mean absolute magnitude in SDSS i-band, and the 
slopes ax of the their composites. 



ri ax z Mi 



1 


-0.84 ±0.04 


2.90 ±0.54 


-27.37 ±0.59 


2 


-0.92 ±0.03 


2.89 ±0.53 


-27.37 ±0.59 


3 


-0.99 ±0.02 


2.87 ±0.51 


-27.36 ±0.60 


4 


-1.02 ±0.05 


2.85 ±0.49 


-27.35 ±0.63 


5 


-1.19 ±0.02 


2.82 ±0.43 


-27.36 ±0.62 


6 


-1.30 ±0.02 


2.84 ± 0.49 


-27.38 ±0.57 


7 


-1.42 ±0.03 


2.82 ±0.44 


-27.41 ±0.63 


8 


-1.42 ±0.02 


2.76 ± 0.40 


-27.43 ±0.57 


9 


-1.49 ±0.05 


2.77 ±0.44 


-27.41 ±0.64 


10 


-1.55 ±0.03 


2.80 ±0.46 


-27.43 ±0.72 


11 


-1.67 ±0.01 


2.75 ±0.45 


-27.39 ±0.67 


12 


-1.83 ±0.04 


2.75 ±0.43 


-27.42 ±0.61 


13 


-1.86 ±0.04 


2.74 ±0.43 


-27.36 ±0.63 


14 


-1.99 ±0.03 


2.72 ±0.42 


-27.30 ±0.61 


15 


-2.00 ±0.02 


2.74 ±0.42 


-27.31 ±0.62 


16 


-2.12 ±0.03 


2.74 ±0.43 


-27.34 ± 0.65 



plotted also diagrams slope ax - absolute magnitude M in 
u (3551 A), g (4686 A), r (6165 A) along with i-band (Fig- 
ure [6|. These magnitudes we calculated using the apparent 
psf magnitudes and reddening in corr esponding bands from 
SDSS DR7, K-correction values from [Richards et all (|2006l ') 
and within the cosmological model mentioned above. As one 
can see from these diagrams there is no dependence between 
the slope and absolute magnitudes in these bands, however 
even the u-band is about 2000 A far from the bump. Hence 
we used also the values of A-' as a characteristic of luminosity 
in this band and plot ax ~ B diagram (Figure [7|, where we 
determined S as a mean flux within 1450 - 1470 A m case 
when the object is located 10 Mpc far from an observer, 
some 'absolute' flux, similarly to definition of absolute mag- 
nitudes. Extending this analogy of absolute magnitudes we 
plotted B in logarithmic scales, but also did not find any 
evidence for dependence of the slope a on luminosity in 
this band. Generally speaking this range is also not the best 
choice, but as it was noticed above this one is the most free 
of obvious emission and absorption features, thus it should 
be the most appropriate to characterize quasar continuum 
luminosity within the bump. 




-3 -2 -1 

slope, a* 



Figure 7. Slope {ax) - 'absolute' flux within 1450 - 1470 A {B) 
diagram. 

3 THE TECHNIQUE 

3.1 General notations 

The data are given in the form of pixels with wavelength 
label Ai, the fl u x de nsity value fj and n oise U j. Following 
iBernardi et all |2003 '): McDonal d~et al.1 (|2006l ') we present 
the observed flux density in i-th pixel within the Lya 
forest region of j-th quasar as 

fl = A^C{\rest){l + (^)(1 + 5^,J + nl (1) 

where Aj is the normalization constant described above 
(Sec. I2.3|l , the wavelength of the absorption Lya feature 
produced by the 'cloud' of intergalactic H I is related to its 
redshift Zi as \i = 1215.67(1 + Zi), C(A) is the mean contin- 
uum level, 8c are deviations of individual continuum from 
the mean one, n is the noise, F is the mean transmission 
for a given redshift and 5f are variance (or fluctuations) of 
transmitted flux. 

Through the fact that different terms of expression ((TJ 
depend on different variables, some of them can be deter- 
mined separately. Namely continuum level C + 5c is a func- 
tion of the rest wavelength Xreat', mean transmission F is a 
function of redshift, 5f depends on a line of sight (presented 
by a given quasar); Ui depends on the instrument properties 
and conditions of observation thus can be considered as a 
function of the observed wavelength Aoi,s; and A-' is a con- 
stant determined in advance. Hence, the mean arithmetic 
composite spectrum for a given subsample with the mean 
redshift z can be presented as: 

f{\rest)^C{\rest)F{z), (2) 
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because the values (Sp^), (5^ J and (n^) are equal to zero 
when averaging over a large n umber of quasar spectr a within 
a given redshift bin (see e. g. iBernardi et al.l l|2003t ) for de- 
tails). Taking into account the similarity of the mean red- 
shifts for all our subsamples we can consider each our com- 
posite spectrum within the Lya-forest region as /(A,est) = 
F'C{\rest), where we assume F' to be the same for each 
subsample. 

3.2 Modelling the spectra 

The mean quasar continuum C{\rest) within the Lya 
forest region can be presented as a sum of the true con- 
tinuum which is usually assumed to have a power-low form 
l|Bernardi et alj [2003 : .Desiacques et al,, .200 7i ) and emission 
lines which are fitted well with Gaussian profile because of 
large velocities of the regions where they are formed. But 
unlike the common treatment we consider the 'pure' contin- 
uum of this wavelength region as a constant, not a power- 
low, due to the fo l lowin g facts. Firstly, it is known (see e.g. 
iDesiacques et af] (|2007l ) for discussion) that extrapolation 
of the power-low 'pure' continuum from the redward part 
of the spectrum yields to underestimated mean transmis- 
sion F in the Lya-forest comparing with the results from 
high-resolution spectra for which continuum fitting within 
Lya-forest region is made directly. This difference is usually 
assigned to large systematic errors, but it is also could be ex- 
plained by the fact that Lya-forest region is located near the 
peak of the Big Blue Bump of quasar spectra, and thus con- 
sideration of the 'pure' continuum within limited wavelength 
band, e. g. redward of 1216 A, as a power-low is just an ap- 
proximation (indeed a tangent to a real curve determining 
UV continuum of quasars, the physical model of which is 
still poorly understood). Secondly, the small width of this 
region together with existence of the emission lines results 
in large uncertainties in parameters of the power-low. That 
is why consideration of such narrow region («150 A wide) as 
a linear (or even a constant in our case) function seems to 
be precise enough. Hence the fitting formula for the product 
of the mean transmission and continuum is the following: 



/(A) = 6 + ^ Ofc exp 



(3) 



Here A is the rest frame wavelength, a^, A^, are the am- 
plitudes, the central wavelengths and FWHM (up to factor 
of \/2) of each fc-th emission feature, and 6 is the constant 
describing continuum. 

Because we are interested more in parameters of the 
lines rather than the global fit of all Lya forest region with 
one function, we split this region onto three parts corre- 
sponding to each emission feature (multiplet) and fit each 
part separately. The fitting procedure was conducted in 
the following two steps: (i) using the IDL Imf it subroutine 
for each multiplet (with central wavelengths around 1070, 
1120 and 1215 A) the best fit model in a form of (|3} was 
found; (ii) the central wavelengths A" obtained at the first 
step were fixed and the best fit values of other parameters 
{6, ak,Wk\ with the la errors from the extremal values of the 
N-dimensional distribution were calculated by the Markov 
Chain Monte Carlo (MCMC) method using CosmoMC pack- 
age as a generic sampler (with the values of {6, afe,iufc} ob- 
tained previously as starting values). The MCMC technique 



has been chosen as it is fast and accurate method of explo- 
ration of high-dimensional parameter spaces. In each case we 
generated 8 chains which have converged to ii — 1 < 0.0075. 

Due to the small values of amplitudes of these lines com- 
paring with the Lya emission line, the central wavelength 
of Lya was determined separately from the red part of the 
third multiplet and then fixed. Note that we used only the 
wavelength range up to 1185 A for this multiplet, thus we ig- 
nored its tiny structure at longer wavelengths and assumed 
the redward part of 1185 A to be only one line with the 
central wavelength about 1215 A. 

Due to the small values the flux dispersion, a'j- , in com- 
posite spectra determined from the covariance matrix to 
estimate errors correctly we introduced additional intrin- 
sic errors, aint, determined such as the total dispersion, 
= CTj is such resulting the minimal x^/d.o.f. to be 

1-0. 



4 THE RESULTS 
4.1 Line identification 

The central wavelengths of all the lines found in com- 
posite spectra are presented in Table 2. Only the lines 
for which the parameters were calculated with the help of 
MCMC method are presented. The arrows stand for regions 
for which the flt failed, 'tent.' means tentative identification 
by eye (which was not included in fit). Despite of variance 
in central wavelengths and the number of lines in different 
spectra we tried to systematize them. In most cases, except 
the multiplet around 1120 A, the lines seem to be the same 
in all spectra, but due to the difference in their FWHM 
in one spectrum two given lines are resolved, while in an- 
other spectrum they are blended and appear to be fitted 
better with one Gaussian. The values in brackets mean such 
dublets fitted by one Gaussian in case of triplets. 

The line identifications were taken from the pa- 
pers on composite UV sp ectra of quasars o b tained with 
HST and FUSE missio ns (|Zheng et all Il997l : iTelfer et al] 
|2002| : IScott et"al] 12004]'), high-resolution s pectra of single 
quasars Ton34 llBinette fc Krongoldl 120081 ) and PHL1811 
(|Leighlv et al.|[2007l ), and al so composite spectr a of quasars 
from optical ground Kas t l|Tvtler et al.l l2004h and SDSS 
jVanden Berk et al.|[200ll ) surveys (see Table g^)). The lines 
which were not identified with known lines from these papers 
are labelled with Xi - Xio. 



4.2 Parameters of lines and continuum 

The obtained values of the line parameters along 
with their la errors from the extremal values of the N- 
dimensional distribution are presented in Tables 4-19, the 
values of constant h for each multiplet are presented in Ta- 
ble (2] The spectra together with the best fits to multiplets 
and the lines found are presented in Figures l9lllOI As one can 
see from Table [2] in most cases the values of continuum pa- 
rameter h within each spectrum vary from one multiplet to 
another within errors and thus the 'true' continuum within 
this wavelength range can be considered as a constant rather 
than having the same power-low form as that within longer 
wavelength range. 
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Table 2. The values of constant b characterizing the 'true' 
continuum for each multiplet. Here n is the number of spec- 
trum. 



uncertainty in estimation of the mean transmission value F, 
is about 20%. 



mult. 1 



mult. 2 



mult. 3 



043+0-020 
q66+° °22 

U.yDD_Q Q4g 
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Figure 8. Ratio of the mean flux, rj^, in each composite spec- 
trum and composite spectrum with ax = —0.84, 770, within the 
wavelength range 1450 - 1470, 1050 - 1100 and 1100 - 1150 A. 



4.3 Normalization constant 

As one can see from increasing of the b values with in- 
creasing of a\ (Sec. 14. 2) ). the first step of the composite spec- 
trum generation, normalization of individual spectra on the 
mean flux within 1450— 1470 A range, introduces additional 
uncertainties into the composite spectra when stacking spec- 
tra with different slopes ax together in one composite. In 
Figure [5] the dependence of the following value 



V = 



{fi)o 



on the number q of the composite spectrum is presented for 
illustration of this fact in Figure 14.21 for wavelength ranges 
1450 - 1470, 1050 - 1100 and 1100 - 1150 A, where the last 
two ranges are located within the Lya forest. Here (fi) is 
the mean fiux within the given wavelength range. As one 
can see from this figure, this deviation of from unity for 
spectra with different slopes, which has to cause additional 



5 DISCUSSION AND CONCLUSIONS 

We presented a new approach for composite spectra 
generation, which lies in stacking only the quasar spectra 
with similar slope ax within the range redward of Lya emis- 
sion line (1215.67 A). The motivation for this is the fol- 
lowing. The quasar spectra in UV range have similar gen- 
eral form: bump-shaped continuum (with maximum around 
1000 — 1300 A), the nature of which is still not clear, with 
broad emission lines similar for all quasars. According to the 
unification scheme of active galactic nuclei this ultraviolet 
bump is formed in accretion disc, hence its shape should re- 
fiect somehow the physical conditions inside the accretion 
disc. Independent on given underlying physics itself one can 
assume that the similarity in shape of one wavelength range 
of the spectrum can imply the similarity in shape of nearby 
ranges. Using own sample of quasar spectra from the 7th 
data release of SDSS we determined the slopes ax of contin- 
uum within the range 1278 — 1480 A considering it to be a 
power-low ~ Ac and generated 16 composites stacking the 
spectra with similar slopes. This technique allowed us to re- 
duce the noise and revealed more tiny features than in case 
of composite spectra made stacking spectra with different 
slopes. 

Considering the wavelength range A ^ 1050 - 1200 A 
known as the Lya-forest region we found more than 14 lines 
within it, three of which were found in previous studies of 
quasar composite spectra from SDSS bv lVanden Berk et al.l 
(2001), some others were found in comp osite spectra from 
space-based HST and FUSE telesc opes (|Scott et al.ll2004l : 
iTelfer et"al]|2002l : lzheng et al.ll 19971 ) or high-resolution spec- 
tra of individual quasars from ground-based telescopes 
iBinette fc Krongoldl l2o"oi : iLeighlv et a"Lll2007l : iTvtler et al] 
[2004). The parameters of the lines were calculated with the 
help of MCMC method. 

Summing up the results we have to note the following. 

• Using the proposed approach helps to reduce the noise 
in composite spectra and study them more carefully. 

• Any dependence of the shape of quasar spectrum within 
the range ~ 1050—1480 A(defined by the slope of continuum 
within 1278 — 1480 A) on luminosity in SDSS u, g, r and i 
bands as far as on luminosity in 1450 — 1470 A band (which 
is the closest to the bump region available for study due to 
absence of emission and absorption lines) was not found. 

• It was shown that the 'true' continuum level within the 
Lya-forest region cannot be considered as a power-low with 
the same slope as in the region redward of Lya emission line. 

The absence of dependence of the UV-bump shape on 
luminosity in the bands mentioned above requires the fur- 
ther study of this region (e. g. more careful study of the 
region redward of Lya emission line) , for understanding the 
physics behind the difference in the UV-bump shape in spec- 
tra of different quasars. It is worth to note that the absolute 
magnitudes in given bands used as a characteristic of lumi- 
nosity include K-correction which is determined within the 
frame of some model of the spectrum general shape. 

The proposed approach can be applied for generation 
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Table 4. Emission lines with their measured rest-frame wavelengths found in individual or composite spectra of quasars within the range 
between hya and Ly/3 emission lines from previous studies. The brackets in case of An and Nil mean that in corresponding paper the 
authors speak about dublet of Arl+Nll. 



line 


Siv 


An 
1066.66 


Nil 
1085.12 


Sim* 
1111.59 


Felll 


Cm* 
1175.67 


Sill 
1194.12 


Sim 


Zhene et al. ( 19971 
Telfer et al. f2002) 
Scott ct al. (2004) 

Binette &; Krongold (2008) 
Leighlv et al. (2007") 
Tytler et^ (2004) 

Vandcn Berk et al. (2001) 


1062+1073 
1067 


tent. 
[1065] 

[1070.95] 
1065.10 


[1084] 
[1084] 
1084.2 


1110.5 


1123 

1123 

1123.13 
1117.26 


1176 

1176 
1175.4 
1175.88 
1175.35 


1195 

1194 
1193.6 


1207 



of new templates for more precise quasar redshift measure- 
ments with the common cross-correlation technique used in 
redshift surveys, more precise theoretical determination of 
K-correction and color-indexes, as far as for determination 
of continuum and mean transmission in Lya forest studies. 
It is also worth to note, that this approach reduces the un- 
certainties in the mean transmission caused by using the 
composite spectra made with common approach, which can 
constitute about 20%. 
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Figure 9. Composite spectra 1—8 with a\ (from top to bottom): 
0.84, -0.92, -0.99, -1.02, -1.19, -1.30, -1.42, -1.42. 



Figure 10. Composite spectra 9—16 with a\ (from top to bot- 
tom): -1.49, -1.55, -1.67, -1.83, -1.86, -1.99, -2.00, -2.12. 
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Table 3. Emission lines found in spectra. Lines in brackets are dublets in case of triplets. Arrows stand for failed fit. 
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Xt 


X8+Cni*+X9 
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1 


tent. 


1058.4+[1071.8] 


1084.9 


1096.3 


1111.8+1125.2 




1161.6+1173.9+1182.1 


1214.4 


2 


tent. 
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1084.6 
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1111.5+1123.3 


tent. 
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Table 5. Parameters of emission lines for spectrum 1 with ax = 



-0.84. 
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Table 6. Parameters of emission 


lines for sped: 


rum 2 with qa 


= -0.92. 
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Table 7. Parameters of emission 


lines for sped: 


rum 3 with a\ 


= -0.99. 
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16.0«:? 



Table 8. Parameters of emission lines for spectrum 4 with a\ 



-1.02. 



Ao, 



1059.5 



1071.0 



1089.3 



1117.0 



1126.2 



1136.5 



1175.2 



1185.3 



1214.7 



096+° °^^ 

U.UyD_o 016 

9 4+0.5 



234+0 045 

U.^oy:_o 032 

q K+i-O 



0.05llHi 



124+0015 
".r^y:_o.oi4 
« q+0.7 

"•'J-0.7 



n 055+0 007 

U.UOO_o 007 
q i:+0.5 



0.084lH}g 
4.7tO;a 



0.091«of, 
5.9i?i 



025+0 055 

U.UZO_o 018 

1 4+2-^ 



1 997+0.157 
l.zz ( _o.o98 

16.lt?i 



Table 9. Parameters of emission lines for spectrum 5 with a\ = 



-1.19. 



Ao,A 1070.7 1087.6 1115.2 1124.4 1136.2 1152.3 1171.9 1214.3 

a 0.22ll°;r8 0.055iroi; 0.043l°:«l« ClOOl^i^v 0.029l«;~? COSSl^fs 0.108t°:«f, 1.122lH«6^ 

w,A S.6t°oi 3.8t°;? 2.2tli S.Otli 2.2«-; 3.2l}:« 7.9l}:J 17.5l«j 

Table 10. Parameters of emission lines for spectrum 6 with ax = —1.30. 

Ao,A 1100.9 1116.9 1125.0 1132.4 1174.2 1213.2 

a u.u^4_Q0(,g u.ii»_(,oog u.U4u_(,oio u.UM_o(,og u.uco_oQi7 i.iiy_oo7o 

W,A O.i_o3 ^-0-0.4 *-'-0.6 J^'-J^_0.6 

Table 11. Parameters of emission lines for spectrum 7 with ax = —1.42. 

Ao,A 1061.6 1071.1 1079.3 1091.0 1108.4 1116.6 1127.0 1143.9 1172.1 1213.2 

a 0.210+o;:?« 0.060+°!;?I 0.137+"°?^ 0.052+o;;« 0.050+oj:^? 0.093+«o?^ 0.117+"°^, 0.042+o«??; 0.084+!!0?5 l.l58+«o^5 



0.019 "•U"U_o (j34 U.iOr_o oig U.UOZ_(, o25 "'UOU.o U.U»0_o pig U.iir_o o23 "■U'±^-0.020 U.UO'i.o oig r--lOO_o o53 

1.5 9/1+0.8 t; Q+1-6 q 1 +2.5 9 S+5-2 S+l-* S S+1-9 q q+3.7 c^/l+1-5 1R9+0-4 

0.9 ^-^-0.9 >J-^-1.0 •J-^-1.3 ^•°-1.2 '^■°-0.S ■J-°-0.8 ■J"^-1.9 >J-^-l.l 10-^_0.4 



Table 12. Parameters of emission lines for spectrum 8 with ax ~ —1.42. 



Ao,A 1060.3 1071.3 1080.8 1098.5 1111.5 1125.4 1141.0 1172.5 1213.8 



a 0.203inf6 0.132troS COSSi^S 0.085lO;g^^ 0.1071^^^ 0.149lHf7 0.082«:«^« 0.112«:«- 



-0.068 



5.2tl:l 3.4«:« 7.7t^-i 3.8lJ:^ 5.0t°-j 6.6t°:« 3.81-? 7.711;? 16.41°;! 



Table 13. Parameters of emission lines for spectrum 9 with ax — 



-1.49. 



Ao,A 



1070.2 



1086.0 



1099.6 



1122.1 



1172.6 



1213.4 



w, 



A 



0.225 
9.9 



+0.007 
0.006 
+0.6 
0.5 



0.034 
1.9 



+0.008 
0.008 
+0.5 
0.5 



0.077 
9.3 



+0.018 
0.016 
+2.1 
1.2 



0.169lHf3 
10.41H 



0.094 

4.8: 



+0.025 
■0.019 
1.2 
0.9 



1 070+0.065 

;+0.5 
-0.5 



17.5: 



Table 14. Parameters of emission lines for spectrum 10 with 



-1.55. 



Ao, A 



1056.3 



1071.7 



1086.8 



1099.6 



1117.3 



1127.4 



1141.7 



1163.7 



1174.3 



1213.1 



a 
w. A 



r, 170+0.030 
U.l ( »_0.020 
q q+0.5 
'^■•^-0.4 



n oo-i +0.029 
« 1 +0-9 



n54+'^ °28 



02^+°-°^^ 

U.U/0_o Qj^4 
9 1 +2.3 
^•^-1.5 



8.otl:» 



0.035lHi? 



024+" °^* 
9 9+3.8 
^■^-1.3 



04fi+° ''i^ 

2-811:? 



nQi+0 012 
u.uyi_Q 

r 7+1.6 



1 1 1 9+0.061 
l.iiZ_Q 

170+0.8 



Table 15. Parameters of emission lines for spectrum 11 with ax 



-1.67. 



Ao, 



1058.1 



1070.4 



1079.9 



1106.7 



1115.7 



1125.1 



1168.5 



1174.9 



1213.8 



r, 990+0.010 
"•^^'-'-0.010 
r, n + 0.9 
' •'5-0.6 



073+'' 0i'5 

A 7+0-7 



0.144lHli 
8.6t°:? 



034+° ''" 
u.uoy:_o 014 

9 4+1.7 



076+° °^2 
u.um_o 010 

n f5+0.7 
■5-0-0.6 



126+''°2^ 

U.iZO_o 013 

6 9+^ '' 



0.069t°:«?° 

4-6t?:^ 



r,P,7+o o22 

U.UO/ _o 026 

2 6+1-'' 

^-0-0.7 



1 91 -1 +0.065 
r-^ii_0.065 
17 0+0.7 
1' --5-0.6 



Table 16. Parameters of emission lines for spectrum 12 with ax = —1.83. 



Ao, A 



1063.0 



1074.1 



1084.0 



1115.5 



1123. i 



1129.0 



1172.9 



1178.5 



2 47+0-041 

W-^^' -0.032 
Q+1.3 
-1.4 



7.3^ 



nS4+'' °^^ 

U.UO4_0.034 
0+0.7 
-0.8 



2.8:* 



0.156 
6.4 



0.041 
0.030 
+ 1.6 
1.3 



r a + 1.5 
0-0-0.8 



040+° °^^ 

U.UW_o.o26 
9 r + 1.1 
^-^-0.9 



0.036tHi^ 

4-ot?:l 



0.076«:«?^ 

4-01?:^ 



0.035«:«^1 



1 r,r,+ 0.064 

-1 ' ' -0.044 

17 3+°-'' 

l'-'5-0.6 



